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Abstract. - We predict the thermodynamic and structural behavior of solutions of long cross- 
linked filaments. We find that at the mean field level, the entropy of self-assembled junctions 
induces an effective attraction between the filaments that can result in a phase separation into 
a connected network, in equilibrium with a dilute phase. A connected network can also be 
Q . formed in a non-thermodynamic transition upon increase of the chain, or cross link density, 

Q_) ' or with decreasing temperature. For rigid rods, at low temperatures, we predict a transition 

from an isotropic network, to anisotropic bundles of rods tightly bound by cross links, that 
is triggered by the interplay between the configurational entropy of the cross-link distribution 
►^ among the rods, and the rotational and translational entropy of the rods. 

(N 

O 

>0 ' Introduction - Branched structures and networks are found in many physical, chemi- 

~^ , cal and biological systems. Physical examples include gels, microemulsions/ wormlike mi- 

^ij ' celles, polymer-colloid dispersions and dipolar fluids [1-5,15,24-26]. Gelation and network 

formation have been studied in a number of experimental [4, 15, 24-27] and theoretical 

C^ ' [7,8,12-14,18,19,28] works and simulations [9,11]. However, the connection between gelation 

and a first order phase separation that commonly accompanies the network formation, remain 

elusive [8-10]. Gels of rigid filaments such as actin cross-Inked with a- actinin [4], or charged 

polyelectrolytes such as DNA, cross-linked with multivalent ions, [21,22] show a qualitatively 

similar thermodynamic behavior, and also have been shown to form compressed bundles of 

^ , filaments, tightly bound by the cross-links. Their origin has been subject of recent theoretical 

activity [20,31]. 

In this paper, we predict, at the mean field level, the phase diagram of a generic system 
of cross-linked equilibrium chains in terms of two independent transitions in the system. The 
first is a continuous, non thermodynamic, topological transition, where an infinite network 
spanning the entire volume of the system is formed. A second, entropically driven, first-order 
thermodynamic phase transition is predicted to occur even in the absence of any specific 
interactions. For most values of the control parameters, this latter transition is a sharp jump 
from a dilute solution of weakly branched chains to a connected dense network. The results, 
shown in Fig^ predict the critical cross-link density and temperature for these transitions. 
In addition, we find that the entropy of the cross links can also drive a transition where 
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the system forms ncm.atically aligned bundles of chains. The bundle transition, as well as the 
phase separation and gelation lines predicted, are consistent with the experiments on actin [4]. 

We study a generic model system that consists of identical monomers that self-assemble 
into chains, far above the polymerization transition, where the chains just start to form [6]. 
The chains, in turn, can branch in the presence of cross links. In order to elucidate the universal 
nature of our predictions, we assume no specific interactions between the monomers, except 
excluded volume, which means that the chains are self-avoiding and mutually avoiding. The 
system is characterized by the following parameters: the monomer density (/), the energies 
of the ends and junctions (branching points) Cg and Cj, measured relative to the energy of a 
bond between two monomers within the chain, the temperature T , the density of cross-links 
c, and the density of the capping molecules responsible for the formation of the ends, p. The 
numbers (per unit volume) of branching points, (pj, and free ends, 0e, are not fixed but are 
predicted as functions of (/>, T, c and p. Here, we focus on the case of sparse junctions and 
ends, (j)j^ 4>e «C 0. The complementary case of very dense junctions is better addressed using 
a different framework and is studied elsewhere [30]. 

This formulation allows us to pinpoint the difference between weak and strong cross- 
linking (i.e., physical gels vs. chemical gels) that correspond to positive and negative junction 
energies e^, respectively. In the case of strong cross links, the density of the cross links must 
be sufficiently low (c ^ (f) so that the total number of junctions formed is small {(t>j < c -^ <j>). 
Of course, in the case of strong cross-linking (ej < 0), the experimental timescale at which 
true, equilibrium thermodynamic behavior is observed, strongly increases as e^ decreases. 

Theoretical Model: We begin by considering the system in the grand-canonical ensemble 
where the network is in equilibrium with a reservoir of ends and junctions of chemical potential 
He and iij respectively. In this formulation the numbers of ends and junctions are not conserved 
but determined by the values of /ig and /ij, respectively; these are calculated below. Both 
junctions and ends can be viewed as "thermal defects" of the system whose "ground state" we 
consider to be an assembly of infinite linear chains. For flexible chains, in the random mixing 
approximation [2], the excluded volume part of the free energy is ■^(j)^. For rigid rods, the 
angular average of the excluded volume interaction between the rods is ■^vn'^, with v ~ dL^, 
n = 4'/L , where d is the rod diameter (taken as unity in our lattice model) and L is the 
average length of the chains [16, 17]. The excluded volume interaction thus scales as ^(j)^ for 
both flexible chains and rigid rods. 

General thermodynamic considerations (as well as rigorous treatments [29]) show that 
in the limit of a low density of non-interacting ends and junctions, every non-conserved, 
thermally generated junction or end lowers the free energy by k^T. At the mean field 
level of analysis, the grand-canonical potential of the system (per unit volume) is given by 
f2(0, /Zj, /ie, Cj, ee) = h<P'^ — <j>e — <t>i whcrc (pe and (fij are the densities of ends and junctions. 
The physics of the problem lies in the dependence of 0e and (j)j on the parameters of the 
system, such as the monomer density (p, temperature T and the chemical potentials jij and 
He- It can be understood via a simple probabilistic argument: two ends are formed whenever 
a bond is broken by thermal fluctuations in the presence of two available end cap molecules. 
To the first order in the number of ends and junctions, the number of bonds is proportional 
to the number of the monomers, (p. Taking into account the free energy of formation of two 
ends, 2(ee — p,e)^ the probability of bond breaking is (/)e^'-''=~'^=''^. In equilibrium, this must 
be equal to the probability of bond formation from the coalescence of two ends. In the 
random mixing approximation, the probability of a collision of two ends is ^gO^ where the 
numerical prefactor ai reflects the microscopic properties of the ends, such as flexibility of 
the chains and effective collision cross-section. We thus flnd 4>e = aie'^''=~^=^/^0^/^. Simi- 
larly, an /-fold junction can form through a collision of / — 2 ends and an internal monomer 
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(or a collision of / ends). The creation of / — 2 ends at the expense of a junction costs 
an energy (/ — 2) (eg — fie) — ^j + l-i-j + Tlnaf, where the coefficient a/ reflects the micro- 
scopic degrees of freedom of the junction, and includes the entropy of configurations of the 
bonds and monomers in the junction. The probability of junction break-up is thus equal to 
</ija/-le(~''^+''j^/'^e(-''~^^(~^=+^'=)/^. Equating this with the probability of collision of / - 2 
ends with an internal monomer, <j)l~'^(j>ai~ , we find for the junction and end densities 

Therefore, in the random mixing approximation, the grand canonical potential of the system 
for the case of sparse junctions and ends {(f>j, (f>e <^ (f>) is : 

f7(</,,M„Me)/r = Id^' - aie(^-^=)/^0i/2 _ afe<^^^-^^y^^f/^ (1) 

Subsequently, we put ai = 1 for convenience. Equation ^ can also be derived rigorously 
[12,28,29] as described in [29], or, in some particular cases, by heuristic methods [11,19,32]. 

It is important to note that at this mean field level it is of no consequence whether in the 
absence of any junctions or ends (zero temperature) the "ground state" consists of infinite 
chains or closed rings. For strong cross-linking, the junctions persist at low temperatures and 
the number of closed linear rings is exponentially small compared with number of branched 
aggregates. Even for weak junctions, the number of closed rings relative to the number 
of chains is smaller by a factor L~^^^ (where L is the chain length). Their influence is of 
importance only at extremely low temperatures and densities, Ref. [6,18]. In addition, the 
role of intra-cluster loops [11] is outside the scope of our model, because the information about 
the long range correlations along a chain is lost in the mean field approximation. 

We now proceed to transform the grand-canonical potential $1(0, /i^) of Eq.(^) and ex- 
press the free energy as a function of the physical cross link and end cap densities, c and 
p, respectively, instead of the chemical potentials /ij and fie- The free energy of a system 
of branched chains in terms of the junction and ends densities (j}j and (j)e, is given by the 
Legendre transform of the free energy Eq. (^ ^(0, (f'jj 4>e) = f^(0, fi-j, Me) + Mj'/'i + f^ed'e'- 

1 /I 

F{(f>,(j)j,(t)e)/T = -<j)'^+(j)j{\n<j)j~l)+<j)e (ln0e - l)+(l)jifj/T-\naf)+(f>eee/T~-(j)j\n(j)--(l)e\n(t) 

(2) 
The first term is the excluded volume of the chains, the next four terms represent the free 
energy of an "ideal gas" of junctions and ends, and the last two terms are the increase in the 
free energy of this "ideal gas" due to the fact that the junctions and the ends are constrained 
to a network. Each junction confines //2 monomers to the same point in space which increases 
the free energy by ^4'j^J■<j> where fi^ ~ Tln0 is the chemical potential of a "gas" of monomers. 
Similarly, last term reflects the reduction in the entropy of a "gas" of ends due to the fact that 
two ends are constrained to lie on the same chain. Because each junction requires one cross- 
link molecule and each end uses up one end cap molecule, the total free energy of the system 

is F*°*((/), 0j, C) = F((/), 0j, 0e) + Faic - <i)j) +Fe{p- (jie) , whcrC Fe{c - (j)j) and Fe {p - (j)e) 

are the free energies of the unbound cross- links and ends, respectively. We take these both to 
be ideal solutions: i^e(^) = Fc (■(/') — T'4}{\nip — 1). The total free energy F*°*((/), (j)j,4>e,c, p) 
is minimized with respect to (j)j and 0e to find the equilibrium end and junction density as 
a function of 4>i ^j-, Se,c, p and T. This predicts that the junction and end densities vary with 
the monomer, cross-link and end-cap densities as 



pe 



=/^0i/V (l + e-^=/^</)i/2^ 0, = ca/e-^^/^0//V (l + a/e'^^/^^^/^) (3) 
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Fig. 1 - (a) The monomer concentration-temperature, I <f), e'^*' /a/ J , plane of the phase diagram for 
strong, four-fold junctions for tj, te < and e4/ee = 1; c = 0.05, p = 0.005. The thick line shows the 
spinodal of the junctions-ends transition. The dashed line is the percolation line to the right of which 
a connected network is formed. The bundles appear to the right of the dotted line; af /q — 1. (b) 
Monomer/crosslink concentration, (0, c), section of the phase diagram for a^e''^''^ — 0.015, |ee/e4| = 
10/3, p = 0.005. The upper line delineates the region of the phase separation. The lower, dashed, 
curve is the percolation line, above which a connected network is formed. Bundles are predicted to 
appear to the right of the dotted line. 



It is important to note that the density of junctions 0^ behaves rather differently for strong 
{tj < 0) and weak (e^ > 0) cross-linking. In the hmit of very strong junctions {ej < 0, |^| 3> 
1), the junction density saturates to a value equal to the cross-link concentration, c. 

This means, that in this limit (?) almost all cross-link molecules are found in the junctions 
{a) the junctions persist at low temperatures. In the opposite limit of weak junctions, (pj ~ 
cafe~^^/'^4'^/'^ <C c which means that (i) most of the cross-link molecules are found in solution 
and not in junctions (ii) at low temperatures, the number of thermally generated junctions 
tends to zero. Also, it follows from Eq.(|3Jl that the condition of sparse junctions, used in the 
derivation of Eq. (Q is always satisfied as long as c <^ 0. 

Substituting the expression of Eq.Q for (pj back into the free energy F, gives: 



F' 



,c) = i02-|-c(lnc-l)+/7(lnp-l)-cln(l+a/e~'^/^(/)^/^)-pln('l-|-e-''=/^0^/^) (4) 



The last two terms are the free energy reduction due to the presence of thermally generated 
ends and junctions. Although both these terms represent negative contributions to the free 
energy, the term associated with the ends is thermodynamically equivalent to a repulsion while 
the term associated with the junctions induces an effective thermodynamic attraction between 

%+<pm+cm-F. 



the monomers. This can be seen by considering the osmotic pressure, II = p^+(t>7n:+c^ 



to which the junctions give a negative contribution, —^cafC 



ajt 



dc 



while the ends give a positive contribution -|-ipe~^=/^0^/^/(l -I- e"'^''/-^^!)^/^). 

Phase Separation and Gelation: If the density of cross link molecules is sufficiently high, 
or the temperature is sufficiently low, the junction induced attraction can be strong enough to 
drive a phase separation between a dense, junction rich, network and a dilute phase consisting 
either of weakly branched chains or a very sparse network. 

For thermodynamic stability, the free energy i^'°*(0,c, p) must be convex with respect 
to all the thermodynamic variables 0, c and p. Technically, this requires that the matrix of 
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second derivatives of i^*°'(0, c, p) is positive definite (i.e., has three positive eigenvalues). This 
condition defines the spinodal surface of the phase separation in the phase space of the system, 
and the critical density of cross-links, Cs{4>, p,T), as a function of the end cap concentration, 
the monomer concentration and the temperature, is given by: 

4 /(/)(4~/)/2 .,/T J ,(i_/)/2 (,^.-c,)/Tx 

In the (c, (j)) plane, the system is stable for cross- link densities c < Cs ((/>), while for c > Cs{(j)) 
the system separates into a dense network with many junctions, in equilibrium with a dilute 
phase. Similarly, in the (T, 0) plane, at a given cross-link density c, the system is unstable 
for a certain range of monomer densities 4> as shown in Fig. 1. The dilute phase can be 
either a very sparse network or consist of disjointed, weakly branched clusters. The dense 
phase is usually a connected network, as discussed below. The thermodynamic behavior is 
rather different for weak and for strong cross-links. In the case of weak cross-links (e^ > 0, 
Cj/T :§> 1) the junction-induced attraction is strong enough to drive the phase separation 
only for three-fold junctions, / = 3. In addition, the phase transition in the case of weak 
junctions is reentrant because the number of thermally generated junctions tends to zero as 
the temperature is lowered as follows from Eq.®. In contrast, in the case of strong junctions 
{cj < 0), the phase separation occurs for junctions of any functionality as follows from Eq.(|$)l: 
the transition is not reentrant because the number of junctions increases with decreasing 
temperature, as follows from Eq.Q . The predicted phase transition is entropic in origin. In 
the dense phase, the entropy of the large number of the formed junctions compensates for the 
loss of the translational entropy of the monomers. 

For sparse junctions, in the mean field approximation [2, 12,29], a connected network 
spanning the entire system is formed when / (/ — 2) (jjj — cj),, . This gelation or percolation 
transition, is a continuous, purely topological transition that has no thermodynamic signa- 
ture. Thus, a connected network spanning the entire system is formed when the cross link 
concentration exceeds the value 



^pcrc 



p/ifif - 2)a/)(/.(i-«/2e(^-'=)/^ (l + e-^^-/^(/.i/2) / (l + aje-'^/^ (bf^') 



This allows one to predict the gelation line in the phase diagram of the system, as shown in 
Fig. 1. 

In vitro experiments on solutions of rigid actin filaments, cross- linked with a-actinin and in 
presence of the capping protein gelsolin, found that the system undergoes a continuous, non- 
thermodynamic transition from an entangled network of branched clusters to a "microgel" [4] 
that results in a sharp increase in viscosity, and can be interpreted in terms of percolation. 
At lower temperatures, it was found that the system undergoes a first-order phase separation 
into a dilute solution of short chains in equilibrium with a dense network that consists mostly 
of bundles of actin filaments, tightly bound by cross-links. These observations are consistent 
with the phase diagram of Fig. 1. 

Bundles: The formation of bundles can be understood in terms of the entropy of rigid chains 
(rods) and cross links that join two rods or four monomers, / = 4. In the bundle phase, the 
rotational and translational entropy of the rods is lower than in the isotropic, network phase. 
However, the entropy of positioning the junctions on the parallel "tracks" formed by the rods 
found in bundles is higher than in the network, where junction positions are confined to the 
intersections between the rods. It turns out that this effect favors the formation of bundles at 
low temperatures. As a first step, we estimate that the transition from an isotropic network 
to bundles takes place when their free energies become equal. Assuming that the energy of 
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the ends is large and negative, so that all the end cap molecules are bound to chain ends, the 
free energy of the isotropic phase (neglecting a term that depends on the number of end caps 
and is the same in both phases) is given by Eq^] 

F,/T ~ i</)2 + c(lnc - 1) - cln (l + afe-'^'^cf'^ + p (in i^p/^^'^) - l) 

We describe the bundle state as a collection of Ni, bundles of m rods of length L. The 
Nj junctions formed by the cross links can be placed between any two adjacent rods. If the 
number of nearest neighbor rods in a bundle is 2q, there are Nb{qm — qm}/"^ /2) available 
one-dimensional "tracks" for the junctions, because the rods at the perimeter of a bundle 
have half as many neighbors as the interior rods; the number of perimeter rods ^ m}'"^ for 
large bundles. The junction entropy is then ~ -^ {qN^Lm [l — 1/ (2to^/^)]) ' . Noting that 
NbLm is the total number of monomers in the system and recalling that the energy of a 
junction is ej, we find for the free energy of the bundle phase 

Fb/T ~ 0, (ln0j - 1) - 0, In (g0 (l - 1/ {^rn^^^))) + ^i^j/T + (c - 0,) (In (c - cp,) - 1) 

where the last term accounts for the translational entropy of the free cross-links in the solu- 
tion. Note that the junction-related part of the bundle state is lower than the corresponding 
part of the free energy of an isotropic network, eq.||2Jl, because the entropic cost of con- 
fining a junction to lie on a track between two parallel rods is —\n(^q(/)(^l — 1/ (2to'^/^))) 
and is lower than the corresponding cost of confinement to an intersection of two rods in a 
three dimensional network, — 21n0. Minimizing the free energy as a function of the junc- 
tion density, (j)j, gives for the junction-related part of the free energy, Fb/T ~ c (Inc — 1) — 
cln (1 + e-^'''^q4> (l - 1/ (2to1/2))). 

In addition to the configurational entropy of cross-link distribution among the filaments, 
the free energy of the bundled states includes the following contributions: («) the translational 

entropy of the bundles ^ I In — -^^^ ^ 1 ) ^-^id (m) the reduction in the rotational entropy of 
the bundles relative to the isotropic phase, p (l — ^) ln47r [16, 17,20] and {in) the excluded 
volume repulsion between the bundles, ^v{(j)/{mL)) where v ~ rn^/'^L'^ [16], because the 
diameter of the bundle of m rods is proportional to m}!'^ . Without (i) and (ii), that balance 
the increase of the configurational cross-link entropy in the bundle phase, the bundles would 
be always preferred thermodynamically. Although the free energy of the bundle state can 
have a minimum at finite values of rn, it turns out that this region is small [32], and for a first 
estimate we assume that the bundle state consists of one macroscopic bundle. In the limit of 
m -^ cxD, FbjT ~ c(lnc— 1) — cln(l -I- e~'^^/-'"0)-|-pln47r. The transition to bundles, therefore, 

takes place when F^ < Fi, that is e'^/"^/a/ < (f) '^'^ ^_/a where Ac = ^(p^ — plnAir, as shown in 
Fig.l. The location of the network to bundle transition line in the phase diagram is sensitive 
to the value of the parameter q/af, where a/ reflects the number of internal configurations of 
a cross-link in the isotropic network relative to parallel bundles; increasing q/af decreases the 
threshold for bundle formation. The parameter a/, as well as the energies eg and ej can be 
estimated experimentally from the number of bound cross-links, and the length distribution of 
chain segments among them. However, when comparing with the experimental results of [4], 
one should bear in mind, that for long chains, the formation of the bundles might be kinetically 
inhibited and may not be experimentally observable. 

* * * 
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